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Abstract 

A new method that relies on computational fluid dy¬ 
namics (CFD) and numerical optimization is used 
to design a transonic business jet wing. The first 
step of this new design method is to develop target 
pressures for a three-dimensional wing design using 
a two-dimensional airfoil optimization code (MSES- 
LINDOP). This airfoil optimization method is fast 
enough to solve a six-point design problem that is 
representative of an entire aircraft mission in a few 
minutes. A full-potential finite element code with a 
solution adaptive cartesian grid (TRANAIR) is used 
to analyze the wing-body-nacelle configuration and 
establish the influence of the fuselage mounted na¬ 
celles on the wing pressures. The blockage in the flow 
caused by these nacelles is approximated in a wing- 
body Euler CFD code (SYN87) with a large bump 
on the aft fuselage. The SYN87 code also solves an 
adjoint set of equations to evaluate the flowfield sen¬ 
sitivities in approximately the same time as that re¬ 
quired to solve the Euler equations for one flowfield. 
These flowfield sensitivities enable three dimensional 
shape optimization in this study with a quasi-Newton 
optimization routine. The objective function used to 
design both the fuselage bump and the wing contours 
was a sum-of-squares of the difference between com¬ 
puted and target wing pressures. Finally, the sur¬ 
face contours are modified slightly with a computer 
aided drawing machine to reduce manufacturing com- 
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plexity. Wind tunnel data from the Boeing Transonic 
Wind Tunnel is in very good agreement with the pres¬ 
sure distributions developed for the 20° swept wing 
considered in this study. This data shows that the de¬ 
sign goals of natural laminar flow at a Mach number 
of 0.75 and minimum wave drag at a Mach number of 
0.80 have been met and provides a validation of the 
design method developed in this study. 


Nomenclature 


Ai 

magnitude of sine bump used to per¬ 
turb airfoil coordinates, used as design 
variable 

c 

wing section chord 

Cl ~ • 

aircraft lift coefficient 

Cp, Cp 

computed and target wing pressure co¬ 
efficients 

eps 

estimate of objective function preci¬ 
sion 

eta 

accuracy parameter for optimization 
line search 

jF 

locations of surface boundary points 

I 

objective or cost function 

Mach , M 

free stream Mach number 

N 

number of pressure coefficients calcu¬ 
lated for a given wing section 

R 

matrix of residuals of flowfield solution 

tol 

convergence tolerance for objective 
function 

V 

vector of design variables 

w 

flowfield quantities, velocity, density, 
etc. 

x/c 

abscissa of normalized wing section co¬ 
ordinate 

X 

locations of flowfield mesh points 

A yjc 

ordinate of perturbation in normalized 
wing section coordinate 

a 

angle of attack 
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^ vector of Lagrange Multipliers 

Subscripts 

i, j counters for spanwise and chordwise 

location of pressure coefficients 
w wing 

Introduction 

Recently, a healthy business jet industry has pursued 
the development of faster aircraft with significant in¬ 
creases in range. 1, 2 ’ 3 Unlike previous business jets 
developed for cruising speeds near a Mach number 
of 0.70, these new aircraft with cruise speeds near 
M = 0.8 require wings with more sweep and carefully 
designed airfoil sections that minimize the onset of 
compressibility drag. Furthermore, these high speed 
wings must be compatible with typical business jet 
configurations that have engines mounted on the fuse¬ 
lage near the wing trailing edge. These requirements 
create an aerodynamic design problem that benefits 
from modern computational fluid dynamics (CFD) 
and aerodynamic shape optimization. 

Several recent developments in aerodynamic sensi¬ 
tivity analysis and the application of CFD to busi¬ 
ness jet design problems make it possible to develop 
a modern transonic wing in a few months. Chan- 
drasekharan, et. al. 4 showed that TRANAIR, 5, 6 a 
full-potential, CFD code with an adaptive cartesian 
grid, could be used to reduce the cruise drag of an 
existing aircraft by tailoring the geometry near the 
wing root and in the nacelle pylon area. Gallman, et. 
al. 7 used a combination of wing-alone optimization 
with R220PT 8, 9 ’ 10 and entire configuration analy¬ 
ses with TRANAIR to design a transonic business jet 
wing for a cruise Mach number of 0.83. Recent Devel¬ 
opments in aerodynamic sensitivity analysis 11 ’ 12, 13 
promise significant improvement in the design pro¬ 
cess. Reference 11 presents the development of a set 
of adjoint equations that enable inexpensive calcula¬ 
tion of flowfield sensitivities. This approach is partic¬ 
ularly appealing because it produces the gradient of 
an aerodynamic objective function with respect to the 
design variables within the same computational time 
and memory as a single flowfield analysis. This is in 
contrast to a finite difference approximation of the 
gradient that requires n 4-1 flowfield analyses, where 
n is the number of design variables. 

This paper presents the application of the SYN87 
design code of Ref. 11 and TRANAIR to the de¬ 
sign of a Mach = 0.80 business jet wing with 20° 
of wing sweep. Figure 1 shows the baseline config¬ 
uration considered in this study. Since SYN87 was 
limited to wing-body configurations, a design process 
was developed that includes the approximation of a 



Fig. 1 Baseline wing-body-nacelle configuration. 


nacelle-pylon configuration with a large bump on the 
aft fuselage. This “bump” approximation was first 
used by Hicks in Ref. 17 to simulate the influence of 
fuselage mounted engines on the wing pressures and 
to avoid the complexity of a multi-block flowfield grid. 
Figure 2 shows an isometric view of a SYN87 wing- 
body-bump model. Full configuration analyses with 
TRANAIR are used to validate the bump approxima¬ 
tion and to design the wing-body fillet. A description 
of the new design process is followed by a discussion of 
the design results, and a comparison with wind tunnel 
data. 

Analysis and Optimization Methods 

The CFD codes used in this study are compatible 
with a fast turn around design process because they 
require simple input files and automate the flowfield 
grid generation. SYN87, used here for wing optimiza¬ 
tion, is based on Jameson’s FL087 14, 15 code, which 
solves the Euler equations using a cell-centered finite 
volume scheme with residual averaging and multigrid 
acceleration. The SYN87 code also includes a wing- 
body grid generator and an efficient calculation of 
the aerodynamic-sensitivity derivatives. TRANAIR 
is a hill-potential, finite-element code with a solution- 
adaptive, cartesian grid. This CFD code requires 
surface-panel geometry as input. TRANAIR was used 
to assess the interference of the nacelle and pylon on 
the wing pressures and to design the surface shape of 
the wing-body fillet. Drela’s MSES and LINDOP 16 
airfoil design design package provided target pressure 
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Fig. 2 Bump approximation of nacelle pylon config¬ 
uration. 

distributions for the three dimensional wing. 

Multipoint Airfoil Optimization 

The interactive airfoil design and optimization frame¬ 
work of MSES/LINDOP enabled the solution of an 
aerodynamic design problem based on six flight con¬ 
ditions. This code solves a coupled Euler/Boudary- 
layer flowfield with a full Newton solution method 
that provides the value of the aerodynamic objective 
function and the flowfield sensitivities or gradients re¬ 
quired for optimization. The full Newton solution 
method provides the gradients with respect to pre- 
defined airfoil design variables for a fraction of the 
cost of one flowfield solution, 16 thereby accomplish¬ 
ing the equivalent of the adjoint solution technique 
used by SYN87. 

Airfoil optimization studies based on objective 
functions defined by a sum of Cd or ML/D for six 
flight conditions identified the compromise between 
the potential for laminar flow at M = 0.75 and the 
requirement for minimum wave drag at M = 0.80. 
The angle of attack and 20 Chebyshev polynomial 
perturbations to the airfoil*coordinates represented 
the design variables that created optimum airfoils 
with a constraint on maximum thickness. Evalua¬ 
tion of the pressure distributions and drag charac¬ 
teristics every 6-8 optimization steps frequently led 
to a slight modification or redefinition of the design 
problem (objective function). Solutions to airfoil op¬ 
timization problems that focused on reducing wave 


drag at M = 0.80, on the development of favorable 
pressure gradients for.laminar flow at M = 0.75, or on 
a weighted combination of these design goals helped 
identify the key trade-offs. A compromise between 
these two aerodynamic goals was then used to define 
the target pressure distributions used in the wing de¬ 
sign process described in the next paragraph. Once 
the three dimensional design process was underway, 
the design-to-target-pressure option in LINDOP was 
used to duplicate the wing pressures in two dimen¬ 
sions. This enabled efficient evaluation of the current 
wing at six design points and guided the development 
of target pressures for 3D wing optimization. 

Wing Design Process 

Figure 3 shows a flow chart of the process used to de¬ 
sign a transonic business jet wing for low wave drag at 
Mach = 0.80 and for natural laminar flow for all Mach 
numbers between 0.75 and 0.80. The MSES/LINDOP 
design tool is used to develop some preliminary air¬ 
foils and to gain a thorough understanding of the six- 
point design problem. A TRANAIR analysis is then 
used to establish the influence of the nacelles on the 
wing and to evaluate the wing pressures. This analy¬ 
sis provides the wing pressure targets for the SYN87 
optimization of a bump on the aft fuselage. The wing 
design problem is now defined with target pressures 
at Mach = 0.80. Airfoil optimization enables the def¬ 
inition of these targets such that the resulting wing 
will have low wave drag at M = 0.80 and the poten¬ 
tial for natural laminar flow at M = 0.75. Typically, 
these targets are also perturbations of the wing pres¬ 
sures obtained in a previous wing optimization. This 
creates a series of design objectives that are small per¬ 
turbations from an existing wing and ensures that the 
targets represent an achievable goal. Finally, SYN87 
is used with a fixed fuselage bump to design a wing 
that minimizes a sum-of-squares measure of the differ¬ 
ence between the computed and target pressure dis¬ 
tributions. At this point it usually becomes neces¬ 
sary to reloft the new wing, and to modify the fairing 
that blends the wing into the body. This lofting step, 
achieved with a computer aided drawing machine, 
reduces manufacturing complexity, ensures that the 
landing gear fit in the wing, and accurately quantifies 
the wing fuel volume. As shown in Fig. 3, the process 
now returns to an analysis of the wing-body-nacelle 
configuration to evaluate the current design. 

Optimization 

The fuselage bump and wing design problems were 
solved using the unconstrained quasi-Newton opti¬ 
mization routine of QNMDIF. 18, 19 The objective 
function for both problems is the sum-of-squares of 


3 



S - s 

' LINDOP-MSES 6 design point 

optimization in 2D. Define 
» 3D target pressure distributions 
v > 


x- : — 1 -- 

' TranAir Analysis 
Mach=0.75,0.80 
^JWing-Body-Nacelle) 


f Syn87 Analysis 
Mach=0.75,0.80 
^(Wing-Body-Bump)^ 


r 




Define new target wing pressures 
at Mach=0.80 for low wave drag 
and natural laminar flow. 

Consider potential for 
laminar flow at Mach=0.75 



/ -\ 

' Syn87 Fuselage Bump Optimization ' 


Minimize difference between 
TranAir & Syn87 wing pressures 


Syn87 Wing Optimization 
Minimize difference between 
computed and target pressures 
using quasi-Newton optimization 


X 


f Reloft wing-body configuration 
with CAD machine to reduce 
manufacturing complexity. 


V 


Fig. 3 Flowchart of design process. 


the difference between computed and target wing 
pressures given by 

a) 

i=l j = l 


where i indicates a specific wing section, j represents 
the chordwise location of the computed and target 
pressure coefficients, and Nj is the number of target 
points at station j. Target wing pressures were defined 
at 14%, 21%, 30%, 47%, 66%, 82%, and 100% semi¬ 
span. The wing surface was constructed by lofting 
linearly between streamwise airfoil sections at 14%, 
21%, 30%. 65%, and 100% semispan. Twenty five 
design variables define the magnitudes of sine wave 
perturbations to the airfoil upper or lower surface at 
each of the wing defining stations. These perturba¬ 
tions to the airfoil coordinates are given by 


25 


i=l 


A y/c =Y, Ai sin 2 ' 7r(x/c) ln ^ 0,5)//In(l / 26) 


( 2 ) 


This type of sine wave was also used to perturb 
the aft-fuselage for the ,: bump :? approximation of the 


nacelle-pylon configuration. Although explicit con¬ 
straints on wing thickness were not included in the ob¬ 
jective function (Eq. 1), a wing fuel volume constraint 
influenced the definition of target pressures through¬ 
out this design study. Each optimization iteration of 
either a wing or bump design problem involves the 
calculation of a desirable search direction using the 
gradient of the objective function and an accumu¬ 
lated approximation of second derivatives, followed 
by a line search through the design space, and a test 
for convergence. Input values for the precision of the 
objective function (eps), an accuracy or convergence 
parameter for each line search (0 < eta < 1), and 
an optimality tolerance (tol = 0.1) for the gradient 
of the objective function control the computational 
tests for convergence in QNMDIF. A fairly accurate 
line search (eta = 0.2) reduced the number of line 
searches required for convergence and improved the 
accumulated approximation of the second derivatives. 
The convergence of the iterative flow solver provided 
an estimate of function precision (eps = 10~ 3 ). For 
the optimization problem considered in this study, 
reductions in the objective function gradient below 
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10” 1 produce imperceptible design improvements or 
variations in the objective function that are less than 
the input function precision. The design process pre¬ 
sented in Fig. 3 does not require convergence to the 
Khun-Tucker conditions for optimality and encour¬ 
ages an assessment of the design improvement after 
every few line searches or design iterations. 

Adjoint Equation Gradients 

SYN87 was used to solve the Euler equations and 
their related adjoint equations on a 193x49x41 mesh 
using the explicit methods of Refs. 14, 15. The ad¬ 
joint equations enable the solution of the objective 
function gradient for an unlimited number of design 
variables. 11 To see how these adjoint equations en¬ 
able the efficient calculation of gradient information, 
consider an objective function given by 

I = I{w,X,T), (3) 

and a set of fiowfield equations represented by 

R(w.X,T) = 0, (4) 

where w represents the fiowfield quantities, X rep¬ 
resents the physical locations of the mesh points 
within the computational domain, and T represents 
the physical location of the boundary. The adjoint 
equation is formulated by combining the gradient in 
the objective function 

rd±\ T — f QL \ T \ du '] 

l dvJ — l dwJ ldv\ 



and the gradient in the fiowfield equations 

qr i _ ra/n ra*.-i 
dv J L dw J L dv J 



After introducing ip, a vector of Lagrange Multipliers, 
Eqs. 5 and 6 are summed 

m T =m T -w T m 

or 

{£} T = {{£} T -w T [|f]}{^} 

+ \m} T -w T m]}m <n 

Since xp in Eq. 7 is arbitrary, it may be chosen such 
that 



is satisfied. The gradient in the objective function 
then becomes 

m T = {{m T -w T m]}m , a . 
+ {{m T -w T m]}m 

where \^\ and {§%} are zero since the surface 
boundary points are considered separately from the 
surface mesh points and the objective function is de¬ 
fined in terms of fiowfield quantities that are evalu¬ 
ated on the surface boundary. All terms except for 
xp are calculated inexpensively using finite differences 
since they only depend on the boundary surface T 
and the mesh X. The computational time required 
to solve Eq. 8 is approximately equal to the time re¬ 
quired to solve the Euler equations for one fiowfield. 
Since a very small amount of computational time is 
then required to solve Eq. 9 for each design variable, 
the total time required for a gradient calculation is 
insensitive to the number of design variables. 

Example SYN87 Wing Optimization 

Approximately 10 loops through the design process 
presented in Fig. 3 produced an acceptable wing de¬ 
sign in this study. However, only two bump opti¬ 
mizations were required to approximate the influence 
of the nacelles on the wing. It is difficult to select 
a wing optimization case from this study that ac¬ 
curately demonstrates the power of SYN87 because 
these wing optimization cases were used to test the 
development of improved boundary conditions for the 
solution, of the adjoint equations. 

An example wing design problem was constructed 
to demonstrate the performance of the SYN87 design 
tool by replacing the defining sections of the final wing 
at 30%, 65%, and 100% semispan with NACA-0012 
sections that have been scaled to the thickness of the 
final wing. The wing sections at 14%, and 21% re¬ 
main unchanged. The design objective was to recover 
the final wing design by specifying the pressure distri¬ 
bution for this wing over the entire span. Fifty design 
variables (Eq. 2) placed at 30%, 65%, and 100% semi¬ 
span define the wing optimization problem. Figure 4 
shows the initial and target pressure at 48% semi¬ 
span. This span station is halfway between the defin¬ 
ing stations at 30% and 65% semispan and the last to 
converge. Figures 5-8 show the level of convergence 
to the target pressures after 3, 5, 7, and 26 optimiza¬ 
tion iterations. After 3 line searches or optimization 
iterations, the shock location and strength match the 
target, but the pressure recovery section is still quite 
different from the target (see Fig. 5). Figure 6 shows 
good agreement between the design and target pres¬ 
sures after just five optimization steps. Seven opti¬ 
mization steps, shown in Fig. 7, take approximately 
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6 hours on a Cray C90 and represent the level of con¬ 
vergence typical for the design process shown in Fig 9. 
Figure 8 shows that 26 optimization iterations and 16 
hours on a Cray C90 produce a nearly perfect match 
between the pressures and recovers the final wing sec¬ 
tion coordinates precisely. 



Fig. 4 Analysis of example wing and target pressures. 



Fig. 5 Convergence in Cp after 3 optimization steps. 


Results 

The MSES-LINDOP airfoil design code enabled the 
development of target pressures that include the com¬ 
promise necessary for acceptable performance at all 
of the flight conditions listed in Table 1. The design 
goals established for these flight conditions include 

• Natural laminar flow in both upper and lower 
surface wing boundary layers for cruise Mach 
numbers greater than 0.75. 

• A minimum amount of nose down pitching mo¬ 
ment 



Fig. 6 Convergence in Cp after 5 optimization steps. 



Fig. 7 Convergence in Cp after 7 optimization steps. 



Fig. 8 Typical convergence in Cp after 26 optimiza¬ 
tion steps. 
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• A limit on the aft loading of wing section pressure 
distributions to enable manual aileron controls. 

• Wave drag characteristics that lead to efficient 
high speed cruise with an overspeed margin of 
0.03 in Mach number. 


• Wing fuel volume sufficient to meet mission range 
requirements. 

• Wing surface contours that facilitate simple low 
cost manufacturing. 

Fortunately, specific constraints were not necessary to 
limit the aft loading and reduce the nose down pitch¬ 
ing moment. Unlike other inviscid airfoil optimization 
studies that lead to significant aft loading in an at¬ 
tempt to reduce the wave drag, the viscous portion 
of the objective function considered by LINDOP cre¬ 
ated more favorable pressure gradients and an accept¬ 
able amount of aft loading. The quarter chord wing 
sweep of the baseline configuration was limited to 20° 
in an attempt to limit the crossflow in the boundary 
layer and increase the potential for natural laminar 
flow. Unfortunately, finite wing effects, flow accelera¬ 
tion caused by the fuselage, and the effective blockage 
created by the fuselage mounted nacelles were suffi¬ 
cient to produce unacceptable wing pressures after 
the initial two dimensional airfoils were used to con¬ 
struct a wing using simple-sweep theory. Figure 9 
shows pressure distributions at the chord break for 
this initial wing both with and without the influence 
of the nacelles. Note the unacceptable shocks on both 
the upper and lower surface and the significant in¬ 
fluence of the nacelle model on the wing pressures. 
Although the results of Fig. 9 suggest that airfoil op¬ 
timization is of limited value in determining the final 
wing design, MSES/LINDOP was invaluable when it 
came to understanding the trade-off between the po¬ 
tential for laminar flow and minimum wave drag. The 
traditional flat-top supercritical pressure distribution 
with a significant amount of aft loading led to lead¬ 
ing edge pressure peaks and adverse pressure distribu¬ 
tions with limited potential for natural laminar flow at 
Mach numbers below the design point. Throughout 
the wing design process presented in Fig. 3, multi¬ 
point airfoil optimization guided the selection of tar¬ 
get pressures. For example, duplicating the 3D wing 
pressures in 2D made it possible to evaluate the in¬ 
fluence of a modification to the wing pressures at 
M=0.80 on the other five flight conditions. 

The extreme difference in the wing-body and wing- 
body-nacelle pressures shown in Fig. 9 emphasizes the 
need to include the nacelles in the aerodynamic analy¬ 
sis. Nacelle inlet and outlet boundary conditions that 
approximate the cruise flight condition and nacelle in¬ 
terference with the fuselage create a large blockage in 


Mach 

Flight Condition 

0.80 

High Speed Cruise 

0.75 

Medium Speed Cruise 

0.70 

Low Speed Cruise 

0.83 

lg Overspeed 

0.80 

1.3g Maneuver 

0.70 

Climb 


Table 1 Flight conditions for multi-point wing opti¬ 
mization. 


.■;*** Wing-Body, M=0.80, C L = 0.44 

-Wing-Body-Nacelle, M=0.80, Cb=0.21 



Fig. 9 Representative pressure distributions for the 
initial wing at a fixed angle of attack. 
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the flow. This blockage reduces the wing lift, even 
near the tip, at a fixed angle of attack and Mach num¬ 
ber. Figures 10-15 show the influence of the nacelle 
blockage on the final wing design. Figures 10 and 11 
show that the nacelle blockage moves the shock for¬ 
ward approximate 20% chord for the inboard portion 
of the wing and Fig. 12 shows a forward shift of 15% 
in shock location at 65% semispan. These figures also 
show excellent agreement between the pressures cal¬ 
culated using the wing-body-bump model of SYN87 
and the wing-body-nacelle model of TRANAIR at ap¬ 
proximately the same lift coefficient. The bump ap¬ 
proximation was found to be quite robust with re¬ 
spect to changes in the free stream Mach number and 
to changes in the incidence of the wing with respect 
to the body, created during the relofting step shown 
in Fig. 3. Figures 13-15 demonstrate the accuracy 
with which the fuselage bump can approximate the 
nacelles at a Mach number of 0.75, even though it was 
designed to simulate the blockage caused by the na¬ 
celles at Mach=0.80. Figures 10-15 also demonstrate 
the insensitivity of the bump approximation to the 
relofting step of the design process, since the TRAN¬ 
AIR model has 0.5° more wing incidence relative to 
the fuselage than the SYN87 model. 


- Syn87, 

Wing-Body-Bump 

C Lw = 0.22 


Wing-Body-Nacelle 

C Lw = 0.23 

-TranAir, 

Wing-Body 

C Lw = 0.34 



i i | i | i | i | i f 

0.0 0.2 0.4 0.6 0.8 1.0 

x/c 


Fig. 10 Wing pressure distributions at 21% semispan 
for optimized wing at Mach=0.80. 

The design objectives of low wave drag and natural 
laminar flow are satisfied with the pressures shown 
in Figures 10-15. At Mach=0.80, Figs. 10-12 show 
favorable pressure distributions to about 45% chord 
on the upper surface and almost 50% chord on the 
lower surface indicating the potential for natural lam¬ 
inar flow. Local Mach numbers less than 1.3 on the 
. wing upper surface produce weak shock waves at ap¬ 
proximately 50% chord along the entire span of the 


-Syn87, 

Wing-Body-Bump 

C Lw = 0.22 


Wing-Body-Nacelle 

C Lw = 0.23 

-TranAir, 

Wing-Body 

C Lw = 0-34 



x/c 

Fig. 11 Wing pressure distributions at 30% semispan 
for optimized wing at Mach=0.80. 


- Syn87, 

Wing-Body-Bump 

Cl w = 0.22 


Wing-Body-Nacelle 

C Lw = 0.23 

- - • TranAir, 

Wing-Body 

C Lw = 0.34 



x/c 


Fig. 12 Wing pressure distributions at 65% semispan 
for optimized wing at Mach=0.80. 
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-Syn87 C Lw = 0.33 

.TranAir Cl w = 0.36 



x/c 


-Syn87 C Lw = 0.33 

.TranAir C Lw = 0.36 



x/c 

Fig. 15 Wing pressure distributions at 65% semispan 
for optimized wing at Mach=0.75. 


Fig. 13 Wing pressure distributions at 21% semispan 
for optimized wing at Mach=0.75. 



Fig. 14 Wing pressure distributions at 30% semispan 
for optimized wing at Mach=0.75. 


wing. This produces a constant shock sweep angle of 
16.5°, even inboard of the chord break where the trail¬ 
ing edge is swept forward. These figures also show a 
slight reacceleration after the shock which could lead 
to boundary layer separation problems, but which will 
also keep the shock from growing as it moves aft with 
increases in free stream Mach number. Figures 13-15 
show the potential for natural laminar flow to ap¬ 
proximately 30% chord on the upper surface and 50% 
chord on the lower surface at a cruise Mach num¬ 
ber of 0-75. Once again the weak shock is located at 
a constant chordwise location, 35%, along the entire 
semispan. 

A 1/8 scale model was tested in the Boeing Tran¬ 
sonic Wind Tunnel to validate our design process. Al¬ 
though the wing contours remain the same, the aft 
fuselage of this model is different from the fuselage 
used in the wing design process and in the analysis 
presented in Figs. 10-15. The improved aft fuselage 
reduced the shock strength between the fuselage and 
the nacelle which creates a slight reduction in the in¬ 
fluence on the wing pressures. The test results pre¬ 
sented in Figs. 16 and 17 for a Mach number of 0.80 
were measured at an angle of attack of 0.88°, and a 
Reynolds number of 2.5 x 10 6 . Figures 18 and 19 show 
the pressures measures for a Mach number of 0.75 at 
angle of attack of 1.4°, and a Reynolds number of 
2.4 x 10 6 . For most of the test, transition disks were 
placed one inch aft of the wing leading edge (approx. 
15% chord), around the aircraft nose, and near the 
leading edge of all other aircraft components. With¬ 
out transition disks on the wing leading edge, flow 
visualization based on infrared thermography showed 
that a laminar boundary layer persisted to the onset 
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of the shock wave. These results provide an assess¬ 
ment of the accuracy of the TRANAIR analyses and 
a validation of the inviscid design methods used in 
this study. 

Figures 16-19 show a comparison of experimental 
data and TRANAIR analyses performed both with 
and without a strip boundary layer on the wing. 
The excellent agreement between the analysis with 
a boundary layer and the experimental data indicates 
that the effects of aeroelastic deflection of the solid 
steel model or interference from the slotted wind tun¬ 
nel walls are negligible. As expected, the bound¬ 
ary layer moves the shock location forward and re¬ 
duces the lift coefficient at the same angle of attack. 
These figures also show that the full-potential shocks 
in the TRANAIR solutions agree with the experiment 
data and indicate that the design goals of weak shock 
waves and low wave drag have been achieved. Fig¬ 
ures 16-17 also show favorable pressure distributions 
with the potential for natural laminar flow, even at 
a slightly higher angle of attack and lift coefficient 
than the design condition presented in Figs. 10-12. At 
Mach=0.75, Figs. 18 and 19 show a slight reduction 
in the potential for upper surface laminar flow and 
a constant chordwise location of a weak shock wave. 
These results validate the wing optimization methods 
based on SYN87 and the TRANAIR calculations. 



Fig. 17 TRANAIR and wind tunnel pressures at 52% 
semispan for optimized wing at Mach=0.80. 


• o Experimental Data, OC = 0.88 
Inviscid, (X = 0.90 

— Viscous, 10% Trip, OC = 0.90 



Fig. 16 TRANAIR and wind tunnel pressures at 21% 
semispan for optimized wing at Mach=0.80. 



x/c 


Conclusions 


Fig. 18 TRANAIR and wind tunnel pressures at 21% 
semispan for optimized wing at Mach=0.75. 


A transonic wing of a high speed business jet was 
.designed using a new iterative design process that in¬ 
cludes wing-body optimization and wing-body-nacelle 
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Fig. 19 TRANAIR and wind tunnel pressures at 52% 
semispan for optimized wing at Mach=0.75. 


analysis. Full-potential analysis of the wing-body- 
nacelle configuration provided an accurate assessment 
of the influence of the fuselage mounted engines on 
the wing pressures. This aerodynamic interference 
was approximated with the wing-body optimization 
code by designing a bump on the aft fuselage. This 
bump was quite effective at capturing the influence 
of the nacelles on the wing pressure distributions and 
enabled the application of the wing-body code to the 
wing-body-nacelle design problem. The wing-body 
optimization tool, SYN87, solves the Euler equations 
using a cell-centered finite volume scheme with resid¬ 
ual averaging and multigrid acceleration. This design 
tool also includes an automatic wing-body grid gener¬ 
ator and an adjoint-equation solver that produces the 
gradient of the objective function in approximately 
the same computational time as the solution of the 
flowfield equations. This approach to the optimiza¬ 
tion problem led to computational times of approx¬ 
imately 6 hours on a Cray C-90 for 7 optimization 
steps or line searches through a design space defined 
by 150 variables. The final wing designed in this 
study satisfies the requirements for low wave drag at 
Mach = 0.80 and the potential for natural laminar 
flow at Mach numbers between 0.75 and 0.80. Pres¬ 
sure distributions measured in the Boeing Transonic 
Wind Tunnel are in excellent agreement with calcu¬ 
lations for the final wing design. These results pro¬ 
vide some validation of the SYN87 wing optimization 
tool and indicate that three dimensional aerodynamic 
shape optimization is now practical for a fast paced 
aircraft development project. 
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